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1. Introduction and statement of the result 

These notes present a new method, developed in collaboration with Johannes Sjostrand 
and Maciej Zworski, the aim of which is a better understanding of quantum scattering 
systems, in situations where the set of classically trapped trajectories at some energy E > 
is bounded, but can be a complicated fractal set. In particular, we are interested in the 
situations where this trapped set is a "chaotic repeller" hosting a hyperbolic (Axiom A) 
flow. Such a scattering system belongs to the realm of quantum chaos, namely the study 
of wave or quantum systems, the classical limit of which enjoy chaotic properties. This 
type of dynamics occurs for instance in the scattering by 3 or more disks in the Euclidean 
plane [15], but also in scattering by a smooth potential (see fig [p. Chaotic scattering 
systems are physically relevant: for instance, mesoscopic quantum dots are often modelled 
by open chaotic billiards [T7] ; the ionization of atoms or molecules in presence of external 
electric and/or magnetic fields also involves classical chaotic trajectories [2J; Open quantum 
billiards can also be realized in microwave billiard expermiments [30J. 

The method we propose is a quantum version of the Poincare section/Poincare map 
construction used to analyze the classical flow (see JQ). Namely, around some scattering 
energy E > we will construct a quantum transfer operator (or quantum monodromy 
operator), which contains the relevant information of the quantum dynamics at this energy, 
in a much reduced form: this operator has finite rank (which increases in the semiclassical 
limit), it allows to characterize the quantum resonances of the scattering system in the 
vicinity of the energy E. The quantum transfer operator is very similar with the open 
quantum maps studied as toy models for chaotic scattering [U [21]. 

Our main result (Theorem [T]) will be stated in §1.21 In §2] and §3] we sketch the proof of 
this result. We defer the details of the proofs, as well as some applications of the method, 
to a forthcoming publication 



From flows to maps, and back. Let us recall some facts from classical dynamics. In 
the theory of dynamical systems, the study of a flow $* : Y — > Y generated by some vector 
field (or ODE) on a phase space Y (say, a smooth manifold) can obten be facilitated by 
considering a Poincare section of that flow, namely a family £ = {Ej, % = 1, . . . , J} of 
hypersurfaces of Y, which intersect the flow transversely. The successive intersections of 
the flow with X define a first return (or Poincare) map k : X — > X (see fig. [I]). This map, 
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Figure 1. Left: a 3-bump potential, which admits a fractal hyperbolic 
trapped set at intermediate energies [251 Appendix]. Right: schematic rep- 
resentation of a Poincare section. 



defined on a phase space E of codimension 1, conveniently represents the flow on Y. Long 
time properties of k are often easier to analyze than the corresponding properties of the 
flow. One can reconstruct the flow from the knowledge of k together with the return 
time function r : E — > IR + , which measures the time spanned between the intersections p 
and k(p). Below we will explain how transfer operators associated with k can also help to 
compute long time properties of the flow. 

1.0.1. Hamiltonian scattering. The flows we consider are Hamiltonian flows defined on 
the cotangent space T*R n . A Hamiltonian (function) p G C°°(T*M n ) defines a Hamilton 
vector field H p on phase space, which generates the flow $* = exp(tH p ). For our specific 
choice (II. ip . the flow is complete. It preserves the symplectic structure on T*IR n , and 
leaves invariant each energy shell p~ l (E), so it makes sense to study the dynamics on each 
individual shell. A Poincare section E C p~ 1 {E) naturally inherits a symplectic structure, 
which is preserved by the Poincare map k. Hence, the Poincare maps we consider are 
(local) symplectomorphisms on E. 

We will specifically consider Hamiltonians of the form 

(i.i) P ( x , e ) = ^! + ^( X ) ; 

with a potential V G C£°(R n ) (say, supported in a ball -6(0, Rq) C W 1 ). Such a Hamiltonian 
generates a scattering system: for any energy E > 0, particles can come from infinity, scatter 
on the the potential, and be sent back towards infinity. Depending on the shape of V and of 
the energy, some trajectories can also be trapped forever (in the past and/or in the future) 
inside the ball -8(0, Rq). This leads to the definition of the trapped set at energy E: 



(1.2) 



K E = {p £p l (E) : exp(RH p )(p) is bounded}, 
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which is a compact, flow-invariant subset of p~ l (E). The interesting long time dynamics 
takes place in the vicinity of Ke, so the Poincare section E need only represent correctly 
the flow <J>* restricted on Ke, or on some neighbourhood of it. The Poincare map k will 

also be defined in some neighbourhood of the reduced trapped set Te = Ke D E. We will 
give a more precise description of E in §2.11 

1.0.2. Chaotic dynamics and transfer operators. Our Theorem [T] will be relevant to the case 
where the flow on Ke is uniformly hyperbolic (and satisfies Smale's Axiom A). Such a flow 
is, in a sense "maximally chaotic". Hyperbolicity means that at each point p £ Ke the 
tangent space T p p _1 (i?) can be split between the flow direction, an unstable and a stable 
subspaces: 

(1.3) T pV -\E) = RH P © E + (p) © E~{p) , 

where the (un)stable subspaces are defined by the long time properties of the tangent map: 
there exist C, X > such that, for any p £ Ke, 

v g E T (p) <=^> Wd^vW < Ce~ xt , t>0. 

The Poincare map k then inherits the Axiom A property. 

To study the long time properties of such chaotic flow, it has proved convenient to use 
transfer operators associated with k [I]. Let us give an example of such operators. Given 
any weight function / £ C(E, R), one define the transfer operator £j by a weighted push- 
forward on functions tp : E — > R: 

Cf<pQ>)= E e /(p 'W)- 

p':K{p')=p 

Provided Cf is applied to some appropriate functional spac^l, its spectrum can deliver 
relevant information about the long time dynamics of k. For instance, the spectral radius 
r sp {Cf) determines the topological pressure of k associated with the weight /, which provides 
statistical information on the long periodic orbits of k: 

log r sp (f) = V(k, f) d = lim — log V" e^ f . 

T— ¥QO 1 * — ' 

| 7 |<T 

(here f is the sum of values of f(p) along the periodic orbit 7). The topological pressure 
of the flow $*, associated with a weight F £ C(X), can also be computed through transfer 

operators. One defines on E the function f(p) = f^ p ^ F($*(p)), that is the accumulated 
weight from p £ E to its next return n(p), and considers the family {jOf_ ST , s £ R} of 
transfer operators. The following relation then relates the pressures of k and 

a = V(¥, F) V(k, f - st) = ^ r sp (£ f _ ST ) = 1 . 



The functional space can be rather complicated, see e.g. [TJ] for the case of Anosov diffeomorphisms. 
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The decay of correlations for the Axiom A flow $* is encoded in the Ruelle-Pollicott res- 
onances, which are the poles of the Fourier transform of the correlation function [22] . 
Within some strip S C C, these resonances {z{\ can be characterized by using the family 
{Cf^ ZT , z G C} of complex weighted transfer operators: Z{ G S is a resonance iff Cf- Zi T 
has an eigenvalue equal to 1. This property can be written (abusively, because transfer 
operators are usually not trace class) as 

(1.4) z G S is a Ruelle-Pollicott resonance •<=>■ det(l — £f- ZT ) = . 

1.1. A quantum scattering problem. Let us now introduce the quantum dynamics we 
are interested in. The operator 

p = p(h) = + v(x) , v g c c °°(ir) , 

generates the Schrodinger dynamics U{t) = exp(—itP(h)/h) on L 2 (IR n ). P{h) is the h- 
quantization of the classical Hamiltonian (II. ip . so the semiclassical behaviour of the quan- 
tum dynamics will be strongly influenced by the flow exp(tH p ). We focus on the dynamics 
around some positive energy E > 0, so the flow we need to understand is $*[ p _1 (i?). We 
will assume that 

• the flow on p~ 1 (E) has no fixed point: dp[ p -ir E \^ 0. 

• the trapped set K E has topological dimension 1. Equivalently, the reduced trapped 
set Te = Ke H £ is totally disconnected. 

These conditions are satisfied, for example, for a 3-bump potential at intermediate energies 
(see fig. [[])• The second condition was absent in previous studies of such systems [281 120] - 
it is a technical constraint specific to the approach we develop below (as we explain after 
Thm[H the condition required for the method to work is actually weaker). 

We are interested in the long time Schrodinger dynamics near energy E, so it is natural to 
investigate the spectrum of P{h) near E. That operator is self-adjoint on L 2 (IR n ), with do- 
main H 2 (M. n ); but, due to the bounded support of V(x), the spectrum of P{h) is absolutely 
continuous on R + , without any embedded eigenvalue. Nevertheless, the truncated resol- 
vent ip(P(h) — z)~ lr ip (with ip G C^°(M n )), well-defined in the quadrant {Re z >, Imz > 0}, 
can be meromorphically extended across the real axis into {Re z > 0, Imz < 0}. The 
finite rank poles {z{\ in this region are called its resonances (they do not depend on the 
specific cutoff ip). Resonances are often understood as "generalized eigenvalues": they are 
associated with metastable modes Ui(x) which are not square-integrable, but satisfy the 
differential equation P{h)ui = Zi Ui, so that they decay expontially in time, at a rate given 
by | ]mzi\/h. 

One of our objectives is to better understand the distribution of these resonances in the 
/i-neighbourhood of the energy E, that is in disks D(E,Ch) (see fig. [2]). More precisely, 
we want to investigate: 
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Figure 2. Schematic representation of the spectrum of P(h) and its reso- 
nances near the energy E. 

• the number of resonances in D(E, Ch). So far fractal upper bounds have been proven 
[28] . We wish to investigate whether similar lower bounds can be obtained, at least 
for a generic system. 

• the width of the resonance free strip in D(E, Ch). A lower bound for such a strip 
has been expressed in terms of some topological pressure rT5J dH [2D] , but a recent 
result of Petkov-Stoyanov (for obstacle scattering) shows that this lower bound is 
in general not sharp [2~T] . 

1.2. Our result. Our main result is a "quantization" of the Poincare section method 
presented above. 

Theorem 1. Assume that, for some energy E > 0, the trapped set Ke for the flow exp(tH p ) 
is topologically one dimensional, and contains no fixed point. 

Then, for h > small enough, there exists a family of matrices {M(z, h), z G -D(0, Ch)} 
holomorphic w.r.to z, such that the zeros of the function 

(1.5) C(z, h) = det(J - M(z, h)) 

give the resonances of (P(h) — E) in D(0,Ch), with correct multiplicities. 

The matrices M(z, h) have the following structure. There exists a Poincare section 
E = Ll/ =1 Sj and map k : X — > S, an h-Fourier integral operator A4(z,h) : L 2 (R n_1 ) J O 
quantizing k, and a projector Uh of rank r(h) x h~ n+1 , such that 

M{z, h) = U h M{z, h) U h + 0(h N ) . 

The remainder estimate holds in the operator norm on C r ( h \ The exponent N can be 
assumed arbitrary large. 

Remark 1.1. The 1-dimensional condition we impose on Ke is not strictly necessary. 
What one needs is the existence of a Poincare section S intersecting <3>* \ Ke, such that 
<9X fl Ke = 0; in particular, we don't need the flow to be hyperbolic on Ke- Still, Axiom 
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A flows provide the most obvious example for which this condition is satisfied [5]; it holds 
as well for the broken geodesic flow in the scattering by 3 disks satisfying a no-eclipse 
condition |15j . 

This theorem shows that the dynamics generated by the Hamiltonian P(h) near E can 
be "summarized" in the family of quantum transfer operators {M(z,h), z 6 D(0,Ch)}. 
One reason for this terminology is that M(z, h) bears some resemblance with the transfer 
operators £f_ 2r briefly described in §1 .0.21 The equation (11.51) characterizing quantum 
resonances is obviously the quantum analogue of the (generally formal) equation (11.41) 
defining Ruelle-Pollicott resonances. Also, the notation ((z, h) in (jl.5p hints at an analogy, 
or relationship, between this spectral determinant and some form of semiclassical zeta 
function (such functions have been mostly studied in the physics literature, see e.g. [7]). 

The operators M(z, h) have the same semiclassical structure as open quantum maps 
studied in the (mostly physical) literature as toy models of quantum scattering systems. 
For instance, the distribution of resonances and resonant modes has proven to be much 
easier to study numerically for open quantum maps, than for realistic flows [H [2H [191 HE]- 
The novelty here, is that the operators M(z, h) allow to characterize a "physical" resonance 
spectrum. 

1.2.1. Historical remarks. Actually, a similar method has been introduced in the theoret- 
ical physics literature devoted to "quantum chaos". To the author's knowledge, the first 
such construction appeared in Bogomolny's work [3] on multidimensional closed quantum 
systems. In that work, a family of quantum transfer operators T(E) is constructed, which 
are integral operators defined on a hypersurface in configuration space. The eigenvalues 
of the bound Hamiltonian are then obtained, in the semiclassical limit, as roots of the 
equation det(l — T(E)) = 0. This work generated a lot of interest in the quantum chaos 
community. Smilansky and co-workers derived a similar quantization condition for closed 
Euclidean 2-dimensional billiards [9], replacing T(E) by a scattering matrix S(E) associ- 
ated with the dual scattering problem. Bogomolny's method was also extended to study 
quantum scattering situations [T2]. On the other hand, Prosen developed an "exact" (that 
is, not necessarily semiclassical) quantum surface of section method to study certain closed 
Hamiltonian systems [23] . 

In the mathematics literature similar operators appeared in the framework of obstacle 
scattering [131 EH] : the scattering problem was analyzed through integral operators defined 
on the obstacle boundaries, which also have the structure of Fourier integral operators 
associated with the bounce map. More recently, a monodromy operator formalism has been 
introduced in [27] to study the Schrodinger dynamics in the vicinity of a single isolated 
periodic orbit. This approach has then been used to investigate concentration properties of 
eigenmodes in the vicinity of such an orbit [5]. The construction we present below heavily 
borrows from the techniques developed in [27]. It improves them on two aspects: first, our 
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invariant set Ke is more complex than a single periodic orbit. Second, the connection we 
establish between the operators {P{h) — E — z) and M(z, h) is deeper than previously. 

2. Formal construction of the quantum transfer operator 

The proof of Thm [Q proceeds in several steps. It uses many tools of /i-pseudodifferential 
calculus (we will use the notations of [HI CD])- We just recall a few of them: 

• a state u = u(h) £ L 2 is microlocalized in a domain U (s T*~R n iff, for any function 
X £ C c °°(T*M n ) with suppxH U = 0, one has ||x™HU 2 = 0(h°°) \\u\\ L 2. (here 
X 1 " = X w ( x i hD x ) denotes the h-Wejl quantization of x)- 

• two states u, v are said microlocally equivalent in U d T*IR n iff, for any cutoff 
X £ C c °°(f/), one has || X W (« - ^)||^ = 0(h°°). 

• similarly, two operators A, B are said microlocally equivalent in V x U (with U, V (s= 
T*M™) iff, for any cutoffs Xi e C c °°(f/), x 2 e C c °°(^), one has \\x2i A ~ b )Xi\\l^v = 
0(h°°). 

• an operator A is microlocally defined in V x U iff it is microlocally equivalent in 
V x U to some globally defined operator. "Microlocally defined in U" will mean 
"microlocally defined in U X IP' . 

The present section constructs the quantum transfer operators microlocally in a neighbour- 
hood of the trapped set K E , without paying attention to the rest of the phase space. The 
arguments making the construction globally well-defined will be presented in §31 

The microlocal construction being strongly tied to a Poincare section, we start by de- 
scribing the latter in some detail. 

2.1. Description of the Poincare section. According to the assumptions of the theo- 
rem, the trapped set Ke is a compact set of topological dimension unity. It is then possible 
to construct a Poincare section S = U^ 1 Sj C p~ 1 {E) with the following properties: 

• each Ej is a (2n — 2)-dimensional topological disk, transverse to the flow. 

• the maximal diameter of the Ej can be chosen arbitrary small. 

• there exists a time r max > such that, for any p £ K E , the trajectory $*(/?) 
intersects E at some time < t < r max . 

• the boundary <9E = UjEj does not intersect Ke- 

If we restrict ourselves to points in the reduced trapped set Te = Ke H S, the map 
p I—)- p+(p) defines a bicontinuous bijection k : T — > T. 

Each component of the reduced trapped set, 71 = f K E H Ej, splits in two different ways: 

(2.1) % = UjDji, where = {p £ % n(p) £ 7^} 

(2.2) % = UjAij, where ^ = {p £ 77, k _1 (p) £ 7^} 
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Figure 3. Schematic representation of a hyperbolic Poincare map. The light 
blue (resp. pink) regions represent .D31 and A 3 i (resp. D 2 i and A21). The 
(un) stable directions are represented by the dashed horizontal and vertical 
lines. The black squares show a coarse-graining of Te- 

We will denote by J+(i) (resp. J-(i)) the set of indices in the "outflow" (resp. "inflow") of 
7i, that is such that Dji and Aji (resp. Dy and Ay) are not empty. The map k is the union 
of components Ky, which relate bijectively Dy with Ay. Since % lies in the interior of E i; 
the components (resp. Ay) are disconnected from one another. Hence, each can De 
extended to be a bijection Ky : Dy — > Ay, where Dy, Ay are open neighborhoods of Dy and 
Aij, respectively in Ej and Ej. The extended map Ky : Djj — >■ Ay is a symplectomorphism 
(see fig. [3] for a sketch). 

2.2. Microlocal solutions. In this section we show that, for z G D(0,Ch), any solution 
to the equation (P(h) — E — z)u = 0, microlocally near some part of Ke, can be "encoded" 
by a transversal function w G L 2 (M n ~ 1 ), which "lives" on one component Ej of the Poincare 
section. 

Take such a component Ej. From the assumption dp \ p -i<e)¥ l O5 there exists an open 
neighbourhood Vj of Ej, and a set of symplectic coordinates (yi, . . . , y n ; , r]i, . . . , r) n ) on Vj, 
such that 

• the Hamiltonian p(p) — E + rji for any p G Vj 

• the section Ej is locally defined by {yi = r]i = 0}, and the origin y = 77 = 
corresponds to a point in %. 

Equivalently, there exists a neighbourhood (0, 0) G Vj G T*IR n and a symplectomorphism 
Ri : Vj — > Vi, such that p o Ki(y, rf) = E + 771 , etc. 

The change of coordinates ki can be /i-quantized into an /i-Fourier integral operator 
Ui : L 2 (M n ) — > L 2 (M. n ), microlocally defined and unitary near Vj x Vj, such that 

(2.3) U* (P(h) — E)Ui is microlocally equivalent with hD yi in Vj x Vj. 
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This "quantum change of coordinates" allows one to easily characterize, for z G D(0,Ch), 
the microlocal solutions to the equation 

(2.4) (P(h) — E — z)u = microlocally in Vj . 
Indeed, the equation 

(2.5) {hD yi - z)v = 
is obviously solved by 

(2.6) v( Vl ,y') = e iz ^ h w(y'), w G L 2 (R n_1 ) , 

that is by extending some "transversal data" w. We denote this extension by v = JC(z)w. 
Conversely, the solution v can easily be "projected" onto the data w. consider some mono- 
tone x £ C°°(M n ) such that x(u) — for yi < e, x(u) = 1 f° r Hi > € - Then, w can be 
recovered from v through 

w(y')= [ e- i ^/ h d yiX (y)v(y)dy 1 . 

JR 

In a more compact form, we write w = JC(z)*x' v, with x' — T[hD yi , X ]- 

The solutions of (12. Ah are then given by selecting some w G L 2 (M n_1 ) (microlocalized 
near the origin), and take 

(2.7) u = U i K{z)w = K: i {z)w. 

That is, the operator /Q(z) builds a microlocal solution of (12 .4p near Sj, starting from 
"transversal data" w G L 2 (M n_1 ). The latter can be interpreted as a quantum state living 
in the reduced phase space The converse "projection" is given by 

(2.8) w = K(z)* X 'U* u = KLiizfx'i u = R +l (z) u . 

Here x% is the cutoff corresponding to X near the section Ej. To get a consistent definition 
for X ii we m ust assume that it jumps back down to a little further along the flow, 
but the precise position will be irrelevant. Indeed, the commutator j^[P(h), xT( x i hD x )} is 
equal (microlocally near Ke) to the sum of two pseudodifferential operators with disjoint 
wavefront sets. The first one is microlocalized near Ej (in the region where X i jumps from 
to 1), we will denote it by x\ = j-[P(h), x7( x i hD x )]i; the second component "lives" in 
the region where X i decreases from 1 to 0, and will not play any role. 

The same construction can be performed independently near each Ej, j = 1, . . . , J. We 
will call Wj the transversal data associated with the section Ej, and )Cj(z), R + j(z) the 
corresponding operators. 
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2.3. Prom one transversal parametrization the next. The the solution ( 12. 7p is mi- 

crolocalized in Vi, since Hi is only defined micro locally in Vi x V{. However, this solution 
can be extended in a forward cylinder Uo<t<T^Sj by using the propagator e -**(- p -- B - z )M : if 
u is a solution near p £ Sj, then e - tt ( p - E - z )/ h u [ s the extension of this solution near 

This way we can extend u up to the vicinity of the sections E 3 - in the outflow of Ej. This 
extended solution will still be denoted by u = JCi(z)wi. Near Ej, this solution can also be 
parametrized by the "transversal" function Wj = R + jU £ L 2 (M n ~ 1 ). The map Wi H- Wj, 
which amounts to changing the transversal parametrization for a single solution u, defines 
our quantum Poincare map: 

(2.9) Mji{z,h) = R +j {z)JCi{z) . 

This operator is a Fourier integral operator quantizing the Poincare map Kjf, it is microlo- 
cally defined, and microlocally unitary, on Dji x Aji. 

Let us see how the operators A4ji(z) can be used. Assume E + z is a resonance of 
P(h), with z £ D(0,Ch). Then, there exists a metastable state u £ Lf oc , global solution 
to the equation (P — E — z)u = 0. The above procedure associates to this solution J 
parametrizations Wi = R + i(z)u, microlocally defined near Tj. For any i and j £ J+(i), 
these parametrizations satisfy Wj = M.ji(z)wi\ can be written in the compact form 

(2.10) w = M(z)w, 

where w = (wi)i is the column vector of all J local parametrizations, and M.{z) is the 
operator valued matrix (Mij(z)). 

In the next subsections we prove that the converse statement holds as well: the existence 
of a solution of (Ai(z) — Id)w = microlocally near Te implies the existence of a solution 
of (P — E — z)u = microlocally near Ke- To prove this we will set up a formal Grushin 
problem, in which the operator (Ai(z) — I) will appear as the "effective Hamiltonian" for 
the original operator (P — E — z). 



2.4. Grushin problems. A Grushin problem for the family of operators^ {(P — E — z) : 
Hl(R n ) — > L 2 (IR n ), z £ D(0,Ch)} consists in the insertion of that operator inside an 
operator valued matrix 

(2.11) V(z) = {^ {P R+ E { ^ Z) R ~f) ■■ H 2 h (M) x H L 2 (R) x H , 

in a way such that V(z) is invertible (see e.g. [29] or [T0~| Appendix] for a general presen- 
tation of this method). Ideally, the auxiliary space "H is "much smaller" than L 2 or H\ (in 



2 -ff^(R") is the semiclassical Sobolev space of norm ||w||#2 = J |u(£)| 2 (l + |/i^| 2 ) 2 <i£, with u the Fourier 
transform of u. 
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our final version, % will be finite dimensional). The inverse of Viz) is traditionally written 
in the form 



v(z) -!_(E(z) E + (z)\ 



The invertibility of (P — E — z) is then equivalent with that of the operator E \-(z): Schur's 

complement formula shows that 
(2.12) 

j(P-E-z)- 1 = E(z)-E + (z)E_ + ( Z y 1 E_(z) , E_ + ( Z y l = -^R + {z){P-E-z)~ l R_(z) 

so that dimker(P — E — z) = dimker E \.{z). For this reason, E ^(z) is called an effective 

Hamiltonian associated with (P(h) — E — z). It has a smaller rank than P(h), but its 
dependence in the spectral parameter z is nonlinear. 



2.5. Our formal Grushin problem. We will first build our Grushin problem microlo- 
cally near Ke (so we can identify H 2 with L 2 ). Our auxiliary space % will contain lo- 
cal "transversal data" Wi G L 2 (R n_1 ), one for each section S», so we have formally 7i = 
L 2 (R n ~ 1 ) J . The auxiliary operators are then vectors of operators: R+(z) = (R+i, ■ ■ ■ , R+j), 
R-{z) = \R-i, • • • , R-j), which will for now be defined only microlocally: 

• R+i{z) is the "projector" ( 12.8ft of L 2 (IR n ) onto the parametrization Wi living on E,. 

We will rebaptize Xi = f Xi (f° r /orward) the cutoff used in the definition of R + i. 

• on the opposite, R_i(z) takes the data W{ G L 2 (IR n ~ 1 ) to produce a microlocal 
solution, and cuts off this solution by applying the derivative of another cutoff \\: 

(2.13) R_i(z) = tfld(z) . 

The cutoff xl (f° r backward) is similar with x{, and xf is? as before, the component 
of [j-P(h), (Xi) w ] microlocalized near Sj. We require that the jump of x\ occurs 
before that of Xn an d that the whole family {x\, i = 1,..., J} satisfies a local 
resolution of identity near Ke' 

(2.14) = ^ ' i R some neighbourhood of Ke- 

i 

2.5.1. Homogeneous problem. Let us now try to invert the matrix V(z) we have just defined, 
at least microlocally near Ke x FJ^ %. First we consider arbitrary transversal data w = (wi), 
and try to solve (in u G L 2 (M™), u_ G L 2 (R n_1 ) J ) the system 

J 

(2.15) Up - E - z)u + ^2R^(z)u_i = 

i=\ 

(2.16) R +i (z)u = Wi, i = l,...,J. 
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Eq. ( I2.16P suggests that u could be a microlocal solution parametrized by u>j, at least in 
the region where Xj jumps from to 1. Since x\ — 1 m this region, we take the Ansatz 

(2.17) u = j^ix^Uz) w t = J2 E +i( z ) w i ■ 

i=l i 

Injecting this Ansatz in ( 12. 15ft . we obtain 

j . j 

(2.18) £ hp, (xtT])Ci(z) + R -i( z ) u -i = » 

i=i i=i 

which we want to solve in (it-i). Each commutator j-[P, [Xi) w ] is the sum of a component 
xf = x[P, (x*i) w ]i microlocalized near Ej, and of components j-[P, ix b i) w ]j microlocalized 
near Aji C Ej, for each index j G J+(i)- The resolution of identity (I2.14p shows that near we 
have H P xf + H P x h j = 0, the quantum version of which reads |[P, {x\) w ]j + ^[-P> (Xj)™]j = 
microlocally near Aj*. As a result d 2 . 1 8 [) can be rewritten as 
j j j 

Near each E^-, j G we have fCi(z)wi = K,j{z)M.ji{z)Wi. For each i we can group 

together the terms localized near Ej, and get: 

R-i(z)wi- 2^ R~i{ z )M ij (z)w j + R_i(z)u^i = Q. 
This leads to the microlocal solution 

(2.19) u^i = -Wi + ^2 -Mij(z) wj = ^2 E~+ij{z)wj . 

We have thus solved the system (12.1 5112. 16p microlocally near Ke x Yl i %, and provided 

explicit expressions for the operators E + (z) and E h (z) — Ai(z) — Id, microlocally near 

the trapped set. 

2.5.2. Nonhomogeneous problem. To complete the microlocal inversion of V(z), we now 
take v G L 2 (IR n ) microlocalized near Ke, and try to solve (in u, w_, microlocally near Ke) 

j 

(2.20) hP-E-z)u + ^2K.i(z)u-i=v. 

i=i 

Let us first assume that v is microlocalized inside the region {x\{p) = l } f° r some index 
i. We then take the truncated parametrix E(z) = J Q e~ lt ^ p ~ E ~ z ^ h dt, with T large enough 
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so that e tTP / h v is microlocalized beyond suppx^, and define the Ansatz u = (xt) w E(z) v. 
The latter satisfies 

(2.21) l -(P - E - z)u = v + I[P, (xtTMz) v dt 

(2.22) =v + tfE(z)v- X-'E(z)v. 

(we have used the splitting of the commutator explained above). Now, provided T is not 
too large, the state E(z)v is microlocalized away from Ej so that xfE(z)v = 0(h°°) On 
the opposite, for each j e J+{i) that state is a microlocal solution of (P — E — z)u — near 
Aji, which can then be written as K,j(z)u_j with u_j = R + jE(z)v. The above equality 
becomes 

— (P — E — z) u = v — R^jU^j , 

and solves f)2.20p microlocally. 

If v is microlocalized near Ej, we cutoff E(z)v by (X^eJ-fi) X^ + Xi)""; which is equivalent 

to identity near 71, and take as above U-j = R + jE(z)v, j G J+(i)- 

We have now fully inverted V(z) microlocally near Ke x J"^ 71, and the norm of the inverse 
can be shown to be of order unity. The effective Hamiltonian reads E_ + (z) = Ai(z) — Id. 
Hence, as anticipated above, the existence of a nontrivial state w satisfying w = M.(z)w 
is equivalent with that of a microlocal solution to (P — E — z)u = near Ke- 

To prove Thm [I], we must define our Grushin problem globally, that is properly define 
the auxiliary space % and the operators R±, in such a way that V(z) is invertible. One 
then says that the Grushin problem is well-posed. 

3. From the formal to the well-posed Grushin problem 

In order to make our Grushin problem well-posed, we will first "deform" the original 
Schrodinger operator P(h) in order to transform its resonances Z\ into bona fide L 2 eigen- 
functions of the deformed operator Pg(h). This deformation is performed through a "com- 
plex scaling" of P(h) far away from the scattering region. The operator {Pe{h) — E) will 
now be elliptic outside a large ball B(0,R). In order to enlarge this zone of ellipticity 
to the complement of a smaller neighbourhood of Ke, we will then modify the Hilbert 
structure of our auxiliary states, using an appropriate escape function G(x,£). After these 
two modifications, we will be able to complete the construction of a well-posed Grushin 
problem, with finite dimensional auxiliary spaces. 

3.1. Complex scaling. Here we use the fact that outside a ball 5(0, Po) 3> supp V, one 
has P(h) = -f ELi &■ In that 

region that operator can be holomorphically exended 
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into P = — ^ Y^k=i acting on functions on C n . For 9 > small, we deform M. n C C n 
into a smooth contour Tg C C n : 

r e n .B C n(o, Rq) = -Br»(o, -R ) > 

T e n C" \ B c « (0, 2R ) = e ie R n n C n \ fi c - (0, 2i2o) . 

We then define the operator Pg(h) acting on u G C^°(Tg), by Pgu = P(u)\r g , where u is an 
almost analytic extension of u. Through the identification Tg 3 x < — > (sin^)" 1 Rex G M n , 
the operator Pg can considered as acting on functions in C^°(M. n ), with the action — e -2 * 61 ^^ 
outside B(0,2Rq). One can then show [26] that the resolvent (Pg — z)" 1 : L 2 — > H\ is 
meromorphic in the region {arg(,z) > —28}. The L 2 spectrum of Pg(h) in that region is 
discrete, independent of 9 and R, and consists of the resonances of the initial operator 
P(h). 

Since Pg(h) = P(h) inside the ball -8(0, Rq) D ttKe, our formal Grushin problem remains 
unchanged if we replace P(h) by Pg(h). 

Below we will take values of 9 of the form 9 ~ C h log(l/7i), C > fixed. 

3.2. Finite dimensional auxiliary spaces. We have built in §2.51 a Grushin problem 
which is invertible microlocally near the trapped set. To make that Grushin problem well- 
posed, we need to make definite choices for the auxiliary spaces, that is for each i — 1, . . . , J 
define a subspace of C L 2 (R n_1 ) containing the transversal data. This subspace should 
contain states microlocalized in some neighbourhood S{ of 71, small enough to lie in the 
domain Uj g j + (j)A,j where k is defined. To construct this subspace explicitly, we may define 
the neighbourhood Si as Si = {<li(p) < 0}, for a well-chosen G satisfying 
liminfp-s.oo qi(p) > 0. The subspace %i can then be defined as the range of the spectral 
projector 

(3.1) Uil R _(qr(y,hD y )). 

According to Weyl's law, for h small enough the space Hi has a finite dimension ~ 
vo\(Si) h~ n+l . 

One can then consider the Grushin problem (12. lip , with P replaced by Pg, the auxiliary 
space H = ® i Hi and the operators 

(3.2) R +i (z) d = f TliK*(z)x b : , R-i = X h :Ki(z)Ui . 

Unfortunately, when trying to solve this new Grushin problem, that is invert V(z), one 
encouters difficulties. Some of them are due to the fact that k does not leave the neigh- 
bourhoods Si invariant (see fig. [3]). As a result, an initial datum u>j G Hi is propagated 
through M.ji(z) into a state M.ji(z)wi which, in general, is not microlocalized in 5^, and 
thus cannot belong to Hj. Brutally applying the projector IL, to Aiji(z)wi produces an 
extra term (1 — Hj)A4jiWi, which is difficult to solve away. Another difficulty arises when 
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trying to solve the unhomogeneous problem (I2.20p for data v microlocalized at some dis- 
tance from Ke- 

3.3. Escape functions and modified norms. These difficulties can be tackled by mod- 
ifying the Hilbert norms on H^iW 1 ) and the auxiliary space L 2 {W l ~ 1 ) J . The new norms 
will be defined in terms of well-chosen escape functions G G C^°(T*M. n ), G 1 G C£°(T*IR n_1 ). 
Using these new norms, the problems mentioned above will disappear, because the states 
microlocalized away from Ke will become easily solvable. 

Let us describe the escape function. For some small 5 > 0, we consider the thickened 
energy shell p~ 1 (E) = IJisi^P -1 ^ + s ) an d trapped set Ke = Ui s |<5 Ke+s- If is shown 
in [2BJ §§4.1,4.2,7.3] and j2UJ §6-1] that, for any small 5q > and large R > 0, and any 
neighbourhoods U C U C V of Ke, one can construct a function Go G C£°(T*M. n ) such 
that 

(3.3) G = on U, H P G > on r£ (0>3H) R n , 

(3.4) H P G >1 onT^rn^yi/), H P G > -5 on T*R n . 

It is convenient!! to slightly modify this function in the neighbourhood of the sets Si C Ej. 
Namely, we consider open neighbourhoods Wi <s Wi of Si in T*M. n , and modify Go into a 
function G h such that H V G X = in Wi while H p d > 1 on T* (0)3ii) R n n(p"-^(5')\(yuU i Wi). 

We then set G = Nh\og{l/h)G\, with iV > fixed but arbitrary large. The exponential 
exp(G w (x, hD)/h) is a pseudodifferential operator in some mildly exotic class, bounded 
and of bounded inverse on L 2 , with norms 0(h~ CN ). We call Hq the vector space Hl(W n ) 
equipped with the Hilbert norm 

(3.5) \\u\\h g = f || exp(-G w (x, hD x )/h)u\\ H 2 . 

Similarly, we consider functions G % G C^°(T*lR n_1 ) such that (using the coordinate change 
hi near E;) G l (y',r]') = Go Ki(0, y'; 0, 77') in some neighbourhood of Si, and modify the 
Hilbert norms on the space L 2 (M n ~ 1 ) attached to the section Ej by 

(3.6) = lle-^^W/Vll^ 1 )- 

3.4. How these norms resolve our problems. Let us explain how this change of norm 
helps us. The action of P$(h) on the Hilbert space Hq is equivalent to the action of 

Pe,c{h) == e~° u ' l h Pe(h) e° w ^ h on if^(IR n ), which is a pseudodifferential operator with sym- 
bol 



%a{p) = p{p) - iNh\og{l/h)H p GM + 0{h 2 \og 2 {l/h)) , P eT B 



(o,RV 



3 The role of this modification is to ultimately keep the norms \\R+i(z)\\H G ^y.i, WR-ii^Wn-j^Ha, 
\Mji(z)\\-Hi^r'Hi uniformly bounded 
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Provided we have chosen a dilation angle 6 <C 5 Nhlog(l/h), the properties of G\ show 
that 

(3.7) Vp £ (V U |J Wi), | Re Pe , G (p) - E\ < 5/2 lmp e , G < -6/C , 

i 

This shows that, for any z G D(0,Ch), the symbol (pe,G — E — z) is invertible outside 
V U Ui^ 7 *' with inverse of order (Mog/i" 1 ) -1 . Hence, for any v G L 2 microlocalized 
outside V U U;^' the equation (Pg.G — E — z)u = v can be solved up to 0(h°°), with 
a solution u microlocalized outside V U Wj. This remark basically tackles the second 
problem mentioned at the end of §3.21 

The first problem (the fact that Kji(Si) is not contained in Sj) is also resolved through 
this change of norms. Indeed, the escape function G\ can be chosen such that it uniformly 
increases (say, by some 2C > 0) along all trajectories of the form p G Si t-> Kji(p) G Ej\ G 

Sj, so that 6 e -G( P )/h — ^ h . This implies that, for any state microlocalized near such 
a point p, one gets (taking the definition (12 M for R+j) 

(3-8) WR+^JQ^WiWH^ = 0(h NC ) \\w t \\ Hcz . 

We then need to modify the finite rank projectors (13. ip defining our auxiliary states, such 
as to make them orthogonal w.r.to the new norms (otherwise ||nj||^ Gi _^H Gi could be very 
large). This modification only amounts to adding a subprincipal (complex valued) term to 
the function q iy such that q™ becomes selfadjoint on H G i, and the spectral projector (13.11) 

orthogonal. The space %i == UiH G i is still made of states microlocalized in Si, and has 
dimension ~ vol(Si)h~ n+1 . Our operators R±i will be defined by (13.21) . 

Let us reconsider the homogeneous problem §2.5.1l with data Wi G Hi in our new Grushin 
problem. For j G J+(i), the state J\4ji(z)iUi does not a priori belong to Hj. However,, 
the estimate (I3.8P shows that the component of Aiji(z) to, microlocalized outside Sj has an 
H G j-norm of order 0(h NC ). As a result, defining the finite rank operators 

(3.9) Mji(z) = Tlj Mji(z) : Hi Uj , 

we find that 

u-i = -Wi + ^ Mij(z)wj G Hi 

provides a solution to the homogeneous problem, up to an error 0(h NC ){^2 i ||u>j||%.). 

The nonhomogeneous problem (I2.20p can be solved as well, up to a comparable error 
(see [IB] for details). 

To summarize, our globally defined Grushin problem has an approximate inverse £(z): 

V{z)£{z) = I + TZ(z), \\n(z)\\ L 2 xH ^ xH = 0{h NC ) , 

where we insist on the fact that N can be chosen arbitrary large (it comes from the factor 
in front of the escape function G). Hence, for h small enough this operator has the exact 
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inverse £(z) = £(z)(I+lZ(z)) 1 = £(z) + L 2 xn ^ H 2 xn (h NC ). In particular, the lower-right 
entry of £ (z) (that is, the exact effective Hamiltonian) reads 

E„ + (z) = I - M(z) + H -> H (h NC ) , 
where M(z) is the matrix composed of the finite dimensional operators (13. 9p . 

As explained in §2.41 this exact inversion implies that the eigenvalues {z{\ of (Pq — E) 
in D(0,Ch) coincide (with multiplicities) with the zeros of det(-EL + (z)). □ 
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